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Abstract 

We present results of statistical analysis of solar wind turbulence using an approach based on the 
theory of Markov processes. It is shown that the Chapman-Kolmogorov equation is approximately 
satisfied for the turbulent cascade. We evaluate the first two Kramers-Moyal coefficients from 
experimental data and show that the solution of the resulting Fokker-Planck equation agrees well 
with experimental probability distributions. Our results suggest the presence of a local transfer 
mechanism for magnetic field fluctuations in solar wind turbulence. 
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I. INTRODUCTION 



Irregular dynamics of the solar wind plasma exhibits many similarities to fully developed 
hydrodynamic turbulence. Numerous in situ measurements of temporal variability of pa- 



rameters of thep 
character [l, 2, 3, 



asma have shown that their spectral distributions usually have power-law 
4J. Investigations of the fluctuations have also revealed their non-Gaussian 
probability distributions at small scales, which is commonly attributed to intermittency 
phenomenon . In fact, the solar wind provides a unique laboratory for studying 



high-Reynolds-number magnetohydrodynamic turbulence (see, e.g., Refs. [2|, [9j for review). 

One of the main problems in the studies of incompressible hydrodynamical turbulence 
is explaining the statistics of velocity fluctuations on different length scales. In magneto- 
hydrodynamic turbulence this problem concerns in general also magnetic field and density 
fluctuations. Conventionally, in investigations of a turbulent cascade, statistical properties 
of fluctuations 5u T (t) — u{t + r) — u(t) of a physical quantity u(t) are examined, where r 
is temporal (or spatial) scale. The fluctuations are studied by examining their probability 
distribution functions (p.d.f.) P(5u T (t)) or n-order moments (5u T (t) n ) of the distributions, 
called also structure functions. Often, if the root-mean-square of velocity fluctuations is 
small as compared to the mean velocity of the flow, one can use the Taylor hypothesis, in- 
terpreting the temporal variation 5u T at a fixed position as a spatial variation 5ui, where / is 
a spatial scale corresponding to the temporal scale r. In an intermittent turbulent cascade, 
the p.d.f. of the fluctuations is non-Gaussian at small scales. When we go to larger scales, 
the shape of the p.d.f. changes, and finally there is a scale t g , such that for t > t g the p.d.f. 



is close to Gaussian distribution 
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111- 



A number of models for the scaling exponents and scaling of the probability distributions 



of the fluctuations have been proposed. Many papers have been a 



tal verification of the proposed models (see, e.g., Refs. 



9. 
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so devoted to experimen- 



ll| for review). Recently, a 



great deal of attention has been devoted to investigations of the fluctuations in hydrody- 



namic 
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turbulence from the point of view of the Markov processes theory (see, e.g., Refs. 



13 



14 



16 



171|). In particular, results of the verification of the validity of the 



Chapman- Kolmogorov equation as well as estimations of the Kramers-Moyal coefficients 
from experimental data suggest that the Markov processes approach may be appropriate to 



the description and modeling of the turbulent cascade 13j, ll4j, ll7(]. The estimations of the 



2 



Kramers-Moyal coefficients allow to determine the form of the Fokker-Planck equation gov- 
erning the evolution of the probability distribution with scale for the fluctuations. A model 
based on a Fokker-Planck equation has been recently proposed for solar wind turbulence, 
but for fluctuations of quantities that exhibit self-similar scaling [18]. In the present paper, 
the Markov processes approach for the first time has been applied to analysis of intermit- 
tent solar wind turbulence. This approach seems to provide a contact point between pure 
statistics and dynamical systems approach to turbulence. 

II. DATA SET 

In the plasma flow expanding from the Sun into the interplanetary space we can distin- 
guish several forms, in particular the slow (< 450 km/s) and fast (> 600 km/s) solar wind 



(see, e.g., Ref. [19| and references therein). At the solar minimum the two forms are usually 
well separated, the fast wind is more homogeneous and incompressible in comparison with 
the slow wind. Our goal here is to study properties of the turbulent cascade, therefore we 
try to exclude effects associated with nonstationary driving and spatial inhomogeneity of 
the turbulence. For this reason, in this paper we have chosen for analysis the fast solar wind 
flowing from non-active high-latitude regions in the solar corona at the solar minimum. This 
data set represents dynamics of fast solar wind free of dynamical interaction with slow wind, 
as possibly the most homogeneous and probably also most stationary case. Therefore effects 
associated with nonstationary driving should be eliminated here to a large extent, and we 
should observe a state possibly closest to freely decaying turbulence, which seems to be 
the most appropriate case to study the turbulent cascade. Since we would like to examine 
fluctuations in a wide range of scales, including small scales, we focus here on magnetic 
field fluctuations, which are available at much better time resolution in comparison with 
measurements of plasma parameters (e.g., bulk velocity or density of the plasma). However, 
we are aware of the importance of detailed analysis of other types of the solar wind, as well 
as other bulk plasma parameters, and we are going to carry out such studies in a future. 

The data set analyzed here consists of about 1.3 x 10 7 measurements of the radial compo- 
nent Br of the magnetic field obtained by the Ulysses spacecraft from 70:1996 to 230:1996 
(day of year:year) at time resolution of one second (see Ref. |20j for the description of the 
experimental setup). The measurements have been obtained at heliospheric latitudes from 



3 



29 to 44 degree and at radial distance from the Sun from 3.5 to 4.2 AU. Small gaps (up 
to three missing points) in the data set have been filled using linear interpolation. Further 
in this paper we consider fluctuations of the radial component of the magnetic field defined 
as b T (t) = Bji(t + r) — B^{t). Presenting results of our analysis, we do not recast the 
fluctuations into the space domain via the Taylor hypothesis, i.e., we consistently use here 
temporal scales. However, since we analyze highly supersonic and super- Alfvenic flow (mean 
velocity U ~ 744 km/s in the reference system moving with the measuring instrument), the 
temporal scales physically should be interpreted rather as spatial scales. Assuming that 
the Taylor hypothesis is satisfied here, one can easily transform the temporal scale r to the 
spatial scale I using the relationship I = Ur [h]]. However, since Ulysses spacecraft provides 
one-point measurements of the magnetic field, in general it is not possible to distinguish 
between temporal and spatial variations in this case. 

In Fig. [1] we show the power spectrum of the radial component of the magnetic field. 
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FIG. 1: Power spectrum of the radial component of the solar wind magnetic field (solid line). The 
dashed line shows the spectrum of the type E(f) oc / _5 ^ 3 for comparison. 



As one can see, the power spectrum has a power-law character with spectral exponent very 
close to -5/3 in the inertial range identified here as stretching approximately from 0.0002 to 
0.075 Hz. 
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III. MARKOV PROCESSES APPROACH 



In the case of a Markov process, by definition the following condition must be satisfied 

P(b Tl , Tx | b T2 , r 2 ; . . . ; b TN , r N ) = P(b Tl , n | b T2 , r 2 ) , ( 1 ) 

thus the iV-point joint p.d.f. P{b T1 . n; b T2 ,T 2 ; . . . ; 6 TJV , rjv) is determined by the product of 
conditional probabilities P(& rj _ 1 , Tj_i|6 Ti , Tj), where 7Vi < r i- For a finite set of experi- 
mental data, the Markov property can be verified by comparison of a conditional p.d.f. 
PE(b T1 ,Ti\b T2 , t 2 ) evaluated directly from data with the p.d.f. computed using the Chapman- 
Kolmogorov equation 

/oo 
P(b Tl , Ti | b T > , r')P(b T > , t'\ b T2 , r 2 )d6 T / , (2) 
-oo 

where T\ < r' < r 2 . Eq. (T5]) is a necessary condition for a stochastic process to be Markovian. 
The Chapman- Kolmogorov equation can be written in a differential form using the so-called 
Kramers-Moyal expansion 

dP(b T ,T\b T0 ,r ) ^ / d x ' 



^\~df) Dik) (br,r)P(b T ,T\b TO ,To). (3) 

U— 1 \ T / 



<9t 

fe=i 

Kramers-Moyal coefficients D^ k \b T ,r) can be evaluated as the limit At — > of the condi- 
tional moments M^(b T , r, Ar), namely 

Z)( fc )(6 T , r) = lim M (k \b T ,T, At) (4) 

A-r^O 

and 

M«(6 T , r, At) = y (6 T , - r'|6 T , t) d6 T ' , (5) 

where At = t — t'. If D^(b T , t) = then according to the Pawula theorem: D^(b T , r) = 
211 ]. In this case, starting from Eq. ([3D we arrive at the Fokker-Planck equation 

dP{b T ,r) ( dD^{b T ,r) d 2 D^(b T1 



for k > 3 



— T- 



V db— + — M — ) P{K ^ (6) 



dr 

which determines the evolution of the probability distribution function of a stochastic process 
generated by the Langevin equation (Ito definition) 



T ^L = D (D (6rj r) + yJ D W(b T ,T) T(t), (7) 
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where T(r) is the delta-correlated Gaussian noise. In comparison with definition used in 



Ref. 



2l| , the Kramers- Moyal coefficients given here are multiplied by r, which is equivalent 



to a logarithmic length scale 171 ] . 

If Eq. §2§ is satisfied, then the transition probability from scale t 2 to T\ can be divided 
into transitions from r 2 to r' and then from r' to T\. Therefore, in the case of a turbulent 
cascade, fulfillment of the Chapman-Kolmogorov equation for all triplets T\ < t' < r 2 in the 
inertial range suggests the presence of a local transfer mechanism in the cascade. 




FIG. 2: Contour plots illustrating verification of the Chapman-Kolmogorov equation for T\ = 750, 
t' = 1000, and r 2 = 1250 seconds. Solid lines represent the conditional p.d.f. PE(b T1 ,Ti\b T2 ,T2) 
evaluated directly from data, whereas dashed lines show the conditional p.d.f. P(6 t1 ,ti|6 T2 ,t 2 ) 
computed using Eq. ([2]). The subsequent isolines correspond to the following levels of the p.d.f.: 
2.0, 0.7, 0.2, 0.07, 0.02 (from the middle of the plot). 



from data and the p.d.f. computed using Eq. (T2]) for n = 750, r' = 1000, and r 2 = 1250 
seconds. One can see that corresponding contour lines for the two probability distributions 
are very close to each other. This indicates that the Chapman-Kolmogorov equation is (at 
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least approximately) satisfied for the range of scales from t% = 750 to t 2 = 1250 seconds. In 
Fig. [3] we show the cuts through the conditional probability distributions for fixed values 



b T =-0.5 [nT] b T =0 [nT] b T =0.5 [nT] 




-1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5 
\ [nT] b X] [nT] b T[ [nT] 



FIG. 3: Verification of the Chapman-Kolmogorov equation (|2j) for t\ = 750, r' = 1000, and 
r 2 = 1250 seconds. Comparisons of cuts through P^{b T1 , T\\b T2 , t 2 ) (points) and P(b T1 ,Ti\b T2 ,T2) 
(lines) from Fig. [2] are shown for fixed values of b T2 , namely b T2 = —0.5 nT (on the left), b T2 = 
nT (in the middle), and b T2 = 0.5 nT (on the right). 



of b T2 . As one can see, points representing cuts through P^(b Tl1 Ti\b T2 , r 2 ) fit well to the 
lines representing cuts through P(b T1 , Ti\b T2 , t 2 ). Repeating such a comparison for different 
triplets Ti, t', t 2 we have checked that Eq. ([2]) is well satisfied in the inertial range (for scales 
from about 50 to 5000 seconds). For larger scales, outside the inertial range, the larger is 
the scale, a worse agreement we observe between experimental p.d.f. and that computed 
using Eq. ([2]). Nevertheless, Eq. ([2]) seems to be fulfilled up to the scale of about 24 hours. 
Therefore, the necessary condition for Markov processes is satisfied here in the entire range 
of scales available for our computations, unlike it is in the case of hydrodynamic turbulence 



as reported in Ref. [17] . where the cascade is not Markovian for small scales, below the 
Taylor length scale. 

We have computed the coefficients (b T , r, At) using the definition of Eq. In Figs. 
IU(a) andH)(b) we present examples of the dependence of the coefficients M^(b T1 r, At) and 
M^ 2 \b T , r, At) on b T for r = 1000 and At = 100 seconds. In Fig. H](c) we show the 
dependence of M^(b T , r, At) on M^(b T ,T, At), which have a more regular and symmet- 
ric form in comparison with the dependence of M^(b T ,r, At) on b T . We propose the 
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FIG. 4: The dependence of the coefficients (a) M^\b T , r, At) and (b) M^ 2 \b T , r, Ar) on 6 r for 
r = 1000 and Ar = 100 seconds. In panel (c) we show the dependence of (b T , r, Ar) on 
M( 1 )(6 r ,r, Ar). 

following approximations: M^(b T , r, Ar) = Al(t, Ar)6 r + A 2 (r, Ar)6^ + A 3 (r, Ar)6^ and 
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M^(b T , r, At) = A A (r, At) + A 5 (t, At)[MW(& t , t, At)] 2 as describing properly the exper- 
imental relationships shown in Figs. S(a) and|4(c), correspondingly. Based on the approx- 
imations, we can fit the parameters Ai(r, At) for a fixed r and changing At, and finally 
compute the limits a* ("7") = liniA T _>o At) for z = 1 ... 5 (e.g., by a linear extrapolation 
toward At = 0) obtaining the following approximations: 

D^(b T ,r) = ai (T)b T + a 2 (T)b 3 T + a 3 (T)b 5 T (8) 

and 

D^(b T .r) = a 4 (r)+a 5 (T)[D^(b T ,r)] 2 . (9) 

Repeating the entire procedure for changing r we can also estimate the dependence of 
the coefficients aj on r. Applying the algorithm, we have obtained the following results 
for the inertial range (r < 5000 seconds) aq = — 3.6t -0,08 , a 2 = 3.5 exp(— O.OOOIt), a 3 = 
— 13.6t~ - 2 , a 4 = 0.00035t 7 , a 5 = 1.2t~ - 3 , and outside the inertial range (t > 5000 seconds) 
a x = -0.5T°- 16 ,a 2 = 2,a 3 = -2.3, a 4 = 0.016t 26 , a 5 = 1.75t-°- 36 . As an illustration, in Fig. 
[5] we show the dependence of the parameter a 4 on r. One can notice a change in the 
dependence for r ~ 5000 seconds, i.e., at the end of inertial range. 




x[s] 

FIG. 5: Dependence of the parameter a 4 on scale r (see Eq. ([9])). 



Parameterization of D^\b T ,r) and D^ 2 \b T ,r) (shown in Eqs. (jHJ) and ([9]) correspond- 
ingly) with experimentally fitted parameters a^r) allows us to solve numerically Eq. (J6]) 
with initial condition taken from parameterization of the experimental p.d.f. at a large scale 
tq, where the probability distribution of fluctuations is approximately Gaussian. Therefore 
we can compute numerically p.d.f. at scales r < tq and compare it to the experimental 
p.d.f., which allows us to verify directly our results. Such a comparison is shown in Fig. [6] 
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FIG. 6: Experimental probability distributions (points) and solution of the Fokker-Planck equation 
(dashed lines) for initial condition (solid line) obtained by approximation of the experimental p.d.f. 
by Gaussian distribution for r = 86400 seconds. We show the comparison of the experimental 
p.d.f. and solution of Eq. for r equal to 86400, 30000, 5000, 1000, 100 seconds (from the top). 
Probability distributions for different scales have been shifted in the vertical direction for clarity 
of presentation. 

for tg = 86400 seconds. As one can see there is a good agreement between experimental 
probability distributions and those computed from the Fokker-Planck equation. 

V. DISCUSSION 

We have shown that the Markov processes approach can be applied to the description of 
the turbulent cascade in fast solar wind turbulence. The Chapman-Kolmogorov equation 
is approximately satisfied in the inertial range, as well as for larger scales up to r = 86400 
seconds. Numerical solution of the Fokker-Planck equation agrees well with experimental 
probability distributions obtained directly from data in the range of r from 100 to 86400 
seconds. Therefore, we can conclude that for intermittent solar wind turbulence, the Markov 
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processes approach can provide a mathematical formalism capable of explaining specific evo- 
lution of the shape of the probability distribution with scale changing from energy containing 
scale down to dissipation scale. Since the formalism describes properly evolution of the prob- 
ability distribution with scale, obviously this should also work for structure functions, which 
are defined as appropriate moments of the probability distributions. Admittedly, direct an- 
alytical derivation of scaling properties of the structure functions can be difficult, but we 
expect that such studies can be done numerically. 

Every Markov process must satisfy the Chapman-Kolmogorov equation (l2lh which ex- 
press the condition that the probability density of transition from scale r 2 to r±, can be 
subdivided into smaller steps, that is transition from scale r 2 to r' , and then from scale r' 
to T\. Therefore, in the case of a turbulent cascade, fulfillment of the Chapman-Kolmogorov 
equation can be interpreted as the presence of a local transfer mechanism in wave vector 
space. The question of locality of the energy transfer is of some interest, e.g., in the studies 
of dynamo mechanism to generate magnetic fields in astrophysical objects, where in helical 
MHD turbulence, nonlocal processes of generation of large-scale fields by small-scale helic- 
ities are studied in details (see, e.g., Sec. 6.2.1 of [HI]). The question is also important for 
modeling MHD flows and numerical simulations, e.g., in large-eddy simulations, where low- 
pass filtering with respect to a cutoff wave number requires some assumptions concerning 
the transfer of energy around the cutoff wavenumber. Local and nonlocal transfer mecha- 
nisms can be distinguished in theoretical studies of turbulence via shell models or numerical 



simulations (see, e.g., [22j, |23|, |24j), but it is very difficult to study the property of turbulence 



using experimental data. The Markov processes approach seems to provide such a method. 
Namely, analyzing a time series from a turbulent flow we should be able to identify the 
character of the dominating transfer mechanism for a given quantity or between different 
quantities, i.e., we should be able to answer the question whether the mechanism is local or 
nonlocal. 

Since our results suggest rather the Markovian character of the turbulent cascade in the 
solar wind, it indicates that the local transfer mechanism dominates in solar wind turbu- 
lence. Therefore dominating transfer of magnetic field fluctuations has similar character 
as in the case of Kolmogorov phenomenology describing turbulence in neutral fluids, where 
according to the picture of Richardson cascade, the energy transfer has local character in the 
wave vector space, i.e., the energy at a scale / is transfered mainly to comparable scales 10]. 
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This result is somewhat surprising, because we analyze here magnetohydrodynamic turbu- 
lence, which is rather magnetic field dominated case. Therefore, according to the classical 
Iroshnikov-Kraichnan picture, due to the Alfven effect, we should expect nonlocal influence 
of large-scale magnetic field on small-scale turbulent eddies, and so also nonlocal interac- 
tions between modes ljj. However, results of recent numerical simulations suggest that 



local transfer mechanisms are present in MHD turbulence 



22 



23 



241 ]. Our paper provides 



experimental results confirming this observation for magnetic-to-magnetic field transfer. 
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